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The second post-Newtonian (2PN) contribution to the 'plus' and 'cross' gravitational wave po- 
larizations associated with gravitational radiation from non-spinning, compact binaries moving in 
elliptic orbits is computed. The computation starts from our earlier results on 2PN generation, 
crucially employs the 2PN accurate generalized quasi-Keplerian parametrization of elliptic orbits 
fSJ , by Damour, Schafer and Wex and provides 2PN accurate expressions modulo the tail terms for 

■ gravitational wave polarizations incorporating effects of eccentricity and periastron precession. 
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I. INTRODUCTION 

Inspiralling compact binaries containing black holes and neutron stars are one of the most promising sources of 
' gravitational radiation for both, almost operational ground based laser interferometric gravitational wave detectors 
>. like LIGO, VIRGO, GEO600 and TAMA300 and the proposed space-based interferometer LISA ||. To obtain an 
CD ■ acceptable signal to noise ratio for detection in the terrestrial detectors, one needs to know a priori the binary's orbital 
. evolution in the inspiral waveform [|j at least upto third post-Newtonian order beyond the (Newtonian) quadrupole 

■ radiation. However, for the measurement of distance and position of the binary, it may be sufficient to know the 
, two independent gravitational wave polarizations and io only 2PN accuracy 0|. Perturbative computation 

I ■ via post-Newtonian (PN) expansions of the binary orbit and gravitational wave (GW) phase are complete to order 
\ beyond the standard quadrupole formula. Extension of the PN perturbative calculations by another two orders, 

^ ' to order w^, is still not complete, because currently used PN techniques |^ leave undetermined a physically crucial 
parameter entering at the level in the gravitational wave flux More recently 0, it has been shown that by 
I . employing several re-summation techniques - to improve the convergence of the PN series - one could make optimal 
use of existing 2PN results to compute GW phasing. Resummed versions of 2PN accurate search templates may be 
. . , just sufficient both for the detection and estimation of parameters of gravitational waves from inspiralling compact 
^ ' binaries of arbitrary mass ratio moving in quasi-circular orbits. For inspiralling non-spinning compact binaries of 
arbitrary mass ratio in quasi- circular orbits, both the 2PN accurate gravitational wave polarizations |^ and the 
associated orbital evolution have been explicitly computed [9|-pl|]. A 2.5PN accurate formula for the orbital phase as 
] a function of time has also been obtained . These expressions are employed by various data analysis packages like 
LAL to search for gravitational waves from inspiralling compact binaries. 

The purpose of the present work is to obtain the 'instantaneous']] 2PN contributions to the two gravitational wave 
polarizations for compact binaries moving in elliptical orbits. On the one hand, these expressions for /i+ and hy 
represent gravitational waves from a binary evolving negligibly under gravitational radiation reaction, incorporating 
precisely upto 2PN order the effects of eccentricity and periastron precession, during that stage of inspiral when the 
orbital parameters are essentially constant over a few orbital revolutions. On the other hand, it is the first (and the 
necessary) step in the direction of obtaining 'ready to use' theoretical templates to search for gravitational waves from 
inspiralling compact binaries moving in quasi- elliptical orbits. The effect of radiation reaction on orbital evolution and 
its consequence on these gravitational waveforms for compact binaries in quasi-elliptical orbits is under investigation 
and will be discussed in the future |Q. 

Galactic binaries, in general, will be in circular orbits by the time they reach the final stage of inspiral. However, 
there exist astrophysical scenarios where compact binaries will have non-negligible eccentricity during the final inspiral 
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^ Following [Q, we term contributions to the GW waveform which depends only on the state of the binary at the retarded 
instant as its 'instantaneous' part. 
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phase. We will next review various such scenarios, some of them speculative, relevant for both ground and space based 
gravitational wave detectors. 

Let us first consider cases that should be important for ground based interferometers. Intermediate mass black hole 
binaries - with total masses in the range 50 Af© < M < (a few) x IQ^Mq - may well be the first sources to be detected 
by LIGO and VIRGO [ p"5Hl^ . Many recent astronomical observations, involving massive black hole candidates point 
to scenarios involving such compact binaries in eccentric orbits. The discovery of numerous bright compact X-ray 
sources with luminosities L > erg / s in several starburst galaxies and rapid time variation of their X-ray fluxes 
implies massive black holes as their central engines. It is suggested that these observations may be explained by the 
merger of globular clusters, containing black holes with M > 10'^ Mq, with its host galaxy However, 10'^ Mq black 
hole present in the center of the globular cluster will have to be created by many coalescences of a > 50 Mq black hole 
with lighter ones and these binaries, in highly eccentric orbits, should be visible to ground based interferometers. This 
scenario may be contrasted with the one suggested in jl^ which also involves compact binaries with high eccentricities. 
However, in this case, black hole binaries, weighing a few solar masses and residing in star clusters, get ejected from 
the cluster by superelastic encounter with other cluster members. These escaping binaries will have short periods 
and high eccentricities before merging. It is worth mentioning that numerical simulations dealing with supermassive 
blackhole formation, performed in the eighties, from a dense cluster of compact stars also indicate creation of short 
period intermediate mass black hole binaries in highly eccentric orbits . 

Recently, there has been studies suggesting that spinning compact binaries may become chaotic [ pO[ . The analysis 
involves numerical evolution of two spinning point masses using 2PN accurate equations of motion. The interesting 
result, observed only for a very restricted portion of the parameter space, is that the outcome of the evolution is highly 
sensitive to initial conditions. It is also observed that binaries whose initial orbits are circular may later become highly 
eccentric. These preliminary results present yet another scenario where eccentricity may become important. 

Many of the potential sources for LISA ||] will be binaries in 'quasi-elliptical' orbits. We list them below, details 
and references to original papers may be found in First, LISA will be sensitive to massive black hole (MBH) 

coalescence involving 10^ to lO^M© black holes, upto 3 Gpc and beyond. It is likely that these binaries will be in 
eccentric orbits during inspiral, as they will be interacting with dense stellar clusters in the galactic nuclei where 
they usually reside. The second candidate involves compact objects orbiting MBH, where compact objects could be 
scattered into very short period eccentric orbits via gravitational deflections by other stars. Finally, LISA will be 
sensitive to thousands of binaries in our galaxy and many of these short period binaries will also be in 'quasi-eccentric' 
orbits. Interestingly, LISA will be highly sensitive to black hole binaries containing primordial black holes of mass 
^ 0.5 A/q. These binaries are one of the speculative candidates for MAssive Compact Halo Objects (MACHOs) [ p2[ . 
It is also shown that the low frequency gravitational waves from black hole MACHO binaries in highly eccentric 

orbits would form a strong stochastic background in the frequency range 10~^ Hz < / < 10^^ Hz, where LISA will be 
most sensitive. 

Finally, we observe that eccentricity will be an important parameter while searching for continuous gravitational 
wave sources in binary systems. Recently, it was shown that searching for gravitational waves from such systems, 
whose locations are exactly known, is computationally feasible p5| . For many such astrophysically interesting systems, 
we note that a post-Newtonian orbital description for generic orbits will be required. 

The computation of the gravitational wave polarizations /i+ and /ix in terms of the orbital phase and frequency 
of the binary was discussed by Lincoln and Will ]26t , using the method of osculating orbital elements from celestial 
mechanics and the 2.5PN accurate Damour-Deruelle equations of motion [^,^. They studied the evolution of 
general orbits and obtained IPN accurate expressions for h+ and for quasi-circular orbits. Later Moreno-Garrido, 
Mediavilla and Buitrago obtained polarization waveforms for binaries in elliptical orbits at Newtonian order with and 
without radiation reaction, studied the effects of orbital parameters and precession on gravitational wave amplitude 
spectrum and implications for data analysis |^,^. Analytic expressions for gravitational wave polarizations and 
far- zone fluxes, for elliptic binaries were obtained to 1.5PN order by Junker and Schafer, and Blanchet and Schafer 
The 2PN accurate gravitational wave polarizations for inspiralling compact binaries moving in quasi- circular 
orbits was given by Blanchet, Iyer, Will and Wiseman For the above calculation they employed the 2PN accurate 
expressions for hjj^ , the transverse traceless part of the radiation field representing the deviation of the metric from 
the flat spacetime and (^), the far-zone energy flux obtained independently using two different formalisms p3|,p|-pT| . 
In the limiting case of a test particle orbiting a Schwarzschild black hole, perturbative calculations are extended to 
very high PN order. For example, in the case of very small mass ratios, polarization waveforms are obtained to 4PN 
order p3. For the case of spinning compact objects in circular orbits, precessional, non-precessional and dissipative 
effects on the gravitational waveform due to spin-orbit and spin-spin interactions were studied extensively [^5|-p8[. 
We note that using the framework we employ here it may be possible to extend results of these papers to compact 
binaries of arbitrary mass ratio moving in elliptical orbits. 

The basic aim of this paper is to obtain the instantaneous 2PN corrections to the 'plus' and 'cross' polarization 



2 



waveforms for compact binaries of arbitrary mass ratio moving in elliptical orbits starting from the corresponding 
2PN contributions to hjj^ • As emphasized in @] , the gravitational wave observations of inspiralling compact 

binaries, is analogous to the high precision radio-wave observations of binary pulsars. The latter makes use of an 
accurate relativistic 'timing formula' based on the solution - in quasi-Keplerian parametrization - to the relativistic 
equation of motion for a compact binary moving in an elliptical orbit In a similar manner, the former demands 
accurate 'phasing', i.e. an accurate mathematical modeling of the continuous time evolution of the gravitational 
waveform. This requires for elliptical binaries, a convenient solution to the 2PN accurate equations of motion. A 
very elegant 2PN accurate generalized quasi-Keplerian parametrization for elliptical orbits has been implemented by 



Damour, Schafer, and Wex |41-44|. This representation is thus the most natural and best suited for our purpose to 
parametrize the dynamical variables that enter the gravitational waveforms. The complete 2PN accurate expressions 
for and consists of the 'instantaneous' contribution computed here supplemented by tail contributions at 1.5PN 
and 2PN orders. The tail computations are not considered here; they must be computed and included in the future. 
The paper is organized as follows: In section^ we present the details of the computation to obt ain ' instantaneous' 



2PN corrections to h+ and hx for inspiralling compact binaries moving in elliptical orbits. Section III deals with the 
influence of the orbital parameters on the waveform. Section ^ comprises our concluding remarks. 

II. THE 2PN GRAVITATIONAL WAVE POLARIZATION STATES 

To compute the two independent gravitational wave polarization states and hy , one needs to choose a convention 
for the direction and orientation of the orbit. We follow the standard convention of choosing a triad of unit vectors 
composed of N, a unit vector along the radial direction to the observer, p, a unit vector along the line of nodes, 
which coincides with y-axis and q, defined by q = N x p (see Fig. 0). The angle between N and the Newtonian 
angular momentum vector which lies along z-axis defines the inclination angle i of the orbit. The orbital phase (f> is 
measured from the positive x-axis in a counter clockwise sense, restricting the values of i from to ^ tt. The two basic 
polarization states and hx are given by 
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h+ ^ 7;{P^P3 - (1^13 ] f^Ij^ ^ (2-la) 



/ix = 2 [P^ Ij + P] ft j Kj > (2-lb) 

where hj^'" is the transverse-traceless (TT) part of the radiation field representing the deviation of the metric from 
the flat spacetime. 



From Eqs. (2.1) it is clear that the explicit computation of 2PN corrections to /i+ and hx requires the fol- 
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lowing: (a) The 2PN corrections to hj^ , generally given in terms of the dynamical variables of the binary, namely 



v'^ , ^^^t't, ni, Wi, N.n, N.v, where r and v are respectively, the relative position and velocity vectors for the two masses 
nil and m2 in the center of mass frame, r = |r|,t) = |v|,n = ^^'f' ~ ^ and m = mi + m2. The unit vector N lies 
along the radial direction to the detector and is given by N = ^, i? being the radial distance to the binary; and (b) 
A 2PN accurate orbital representation for elliptical orbits to parametrize these dynamical variables. 

Before explaining in detail the procedure to compute 2PN contributions to h+ and hx, we will first illustrate that 
computation by presenting in detail the Newtonian computations for /i+ and hx- 

A. The Newtonian GW Polarizations 

At the leading Newtonian order, we have 



iKJn)N = ^^yfem(N) l^V., - — 71,, j , (2.2) 

where Pijkmi^) is the usual transverse traceless projection operator projecting normal to N, Vij = ViVj , riij = riirij; 
and /i is the reduced mass of the binary, given by mi m2/m. Note that the above contribution arises from the mass 
quadrupole moment of the binary. 



There is no need to apply the TT projection in Eq. (|2.2|), and Eq. (2.1) at the leading order gives 
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hx = 
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(2.3a) 



(2.3b) 



The convention we adopted to define the triad of unit vectors implies p = (0,1,0), q = (— cos i, 0, sini), N = 



(sini, 0, cosz ), n — (cos0, sin0, 0), and -v — [f cose 
inputs, Eq. (|2.3|) becomes, 



sin0, f sin0 + r (p cos0, 0), where — d(j)/dt. With these 
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sin 2(j) — 2r r (p cos 20 



(2.4a) 
(2.4b) 



where rj = n/m and C and S are shorthand notations for cosi and sinj. 

When deahng with eUiptical orbits, it is convenient and useful to use a representation to rewrite the dynamical 
variables r, r, (p and (j) in terms of the parameters describing an elliptical orbit. For example, in Newtonian dynamics, 
the Keplerian representation in terms of eccentricity, semi-major axis, eccentric, real and mean anomalies is a conve- 
nient solution to the Newtonian equations of motion for two masses in elliptical orbits. The Keplerian representation 
reads: 



r = a (1 — e cos u) 
n(t — Iq) = I = u ~ e sin u , 
(p- (po^ V, 

where v — 2 tan^ 




tan(-) 



(2.5a) 
(2.5b) 
(2.5c) 

(2.5d) 



where u, Z, v are the eccentric, mean and real anomalies parametrizing the motion and the constants a, e, n, to, (pQ 
represent semi-major axis, eccentricity, mean motion, some initial instant and the orbital phase corresponding to that 
instant respectively. These constants which characterize a given eccentric orbit may be expressed, at the Newtonian 
order, in terms of the conserved energy E and angular momentum per unit reduced mass | J| as 



(2.6a) 
(2.6b) 
(2.6c) 



(-2i?)' 
e^l + 2Eh'^ , 

""^^^^ 

with h = Note that n — where T is the orbital period. 

In the case of circular orbits e = and v — u — I, (p is thus a linearly increasing function of time and r = 0, 
(j) — n — The polarizations are uniquely given by the straightforward substitutions of these simple limiting forms. 
The only residual choice is whether one uses the gauge-dependent variable 7 = or the gauge-independent variable 
X = (vrmi^Gw)^^'^- The situation is more involved in the case of general orbits even at the leading Newtonian order. 
Indeed, if e ^ 0, then v u ^ I, v(u) and u are more complicated functions of I and thus (p is not a simple linearly 
increasing function of time. This is why a straightforward representation of the polarizations in terms of v and u or 
even a more involved one in terms of u only, which may be obtained by explicit elimination of v are inadequate. The 
clue to the correct description follows from the analysis of Damour Bs] for the 2PN accurate equations of motion of a 
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compact binary. Here it was shown that the basic dynamics can be represented as a function of two variablesg A and I 
and be 2tt -periodic in both of them. The GW polarizations will inherit this double periodicity and we shall crucially 
exploit it as follows: We will split (j) into a part A linearly increasing with time and the remaining part denoted by 
W(^) which is a periodic function of I: 



= A + W(0, 



where A and W(^) are given by 



W(0 = {v{u{l)) 



u{l) + e sinM(Z)) . 



(2.7) 



(2.8a) 
(2.8b) 



It is worth emphasizing following [[45[ that though one could consider A to be a function of / given by A = 0o + ^ it is 
more advantageous to consider (f) and consequently the GW polarizations to be independently periodic in both A and 
I. With this observation, as we will see below, it is natural to split the A and W(^) dependence in the polarization and 
consider the A dependence as representing the harmonic time dependence and the W(Z) term as representing the time 
varying amplitude modulation ('nutation'). This clean separation also facilitates a simple and precise treatment of 
the spectral decomposition of the GW polarizations as shown in Section [11 . Finally, we note that this decomposition 



is not just appropriate to discuss effects of eccentricity on the Newtonian waveform and the periastron precession at 
IPN order but powerful enough to analyze all PN effects upto the 2PN order. 

Armed with the above important conceptual input, the computation of hj^ and involves a routine, albeit lengthy 
algebra. It is straightforward to obtain Newtonian expressions for r, r and (p in terms of n, e and u using Eqs. ( p.E| ), 
( p.6D and relations for E and h^, given at Newtonian order by {—2E) = {Gmn )i, {-2Eh'^) = 1 - e^. They are 
given by 



Gm 1 

r = ( — — ) 3 (1 — e cosm) , 

(Gmn) 3 e sinu 

r ~ 

(1 — e sinu) ' 

^ n \/l ~ 

(1 — e sinw)2 



(2.9a) 
(2.9b) 
(2.9c) 



Using Eqs. (2.9) and splitting </) as </> = A + W(/) in Eqs. (2.4), we obtain after some manipulation: 
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Gmrj .Gmn 2/3 

R (1 ~ e cos uY 



(l + G^) 



(e cosw) — e cosm 



2e^-t-2 cos2W + 2e\/l - sinu sin2W 



cos 2 A 



+ (1 + G2) 



2 e ^/l — sin u cos 2 W — ( (e cos u)^ — e cos u 



2 6^ + 2 I sin2W 



Gmrj .Gmn 



sin 2 A + S'^ e cos u{l — e cos u) 
1 



(2.10a) 



■G 



R (1 ~ e cos u)2 

+ ( 2 (e cos u)^ - 2 e cos u - 4 + 4 ) sin 2 W 



— ( 4 e Vl — sin u ) cos 2 W 
cos 2 A 



2 (e cos u) - 2 e cos -u - 4 e + 4 cos 2 W 



4 e \/l — sin u ] sin 2 W 



sin 2 A 



(2.10b) 



We denote by A the variable denoted by m in Ref. 
including here. 



to avoid confusion with the total mass m in most current literature 
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We observe that terms like sin 2 W sin 2 A, c os 2 W cos 2 A, and cos 2 W sin 2 A, sin 2 W cos 2 A are generated respec- 
tively from cos20 and sin20 terms in Eqs. (2^) due to the cj) = W + X split. The circular limit of Eqs. ( 2.10 ), 
obtained by putting e = 0, agrees with the Newtonian terms in Eqs. (2), (3) and (4) of 



B. The 2PN GW Polarizations 



The computation of 2PN corrections to /i+ and ft-x is similar in principle to the Newtonian calculation. However, 
there are subtleties and technical details which will be presented, in some detail, below. 

From the Newtonian calculations, it is easy to note that we require a 2PN accurate orbital representation for 
computing 2PN corrections to h+ and hx- We employ the most Keplerian-like solution to the 2PN accurate equations 
of motion, obtained by Damour, Schafer, and Wex [^2|-Q, given in the usual polar representation associated with the 
Arnowit, Deser and Misner (ADM) coordinates. It is known as the generalized quasi-Keplerian parametrization and 
represents the 2PN motion of a binary containing two compact objects of arbitrary mass ratio, moving in an elliptical 
orbit. The relevant details of the representation is summarized in what follows. 

Let r{t), (pit) be the usual polar coordinates in the plane of relative motion of the two compact objects in the ADM 
gauge. The radial motion r{t) is conveniently parametrized by 

7' = flj. (1 — cosw) , (2.11a) 
f 9 

n(t — to) — I — u — et sin u -\ — j sin v -\ — j (v — u) , (2.11b) 

where u is the 'eccentric anomaly' parametrizing the motion and the constants Gr, e^, et, n and to are some 2PN 
semi-major axis, radial eccentricity, time eccentricity, mean motion, and initial instant respectively. The angular 
motion (f){t) is given by 

l-|-^^w-f-^sin2i;-l-^sin3w, (2.12a) 
where v = 2 tan"! | (^^-^) ' tan(|)| . (2.12b) 

In the above v is some real anomaly, (j)Q,k,e^ are some initial phase, periastron precession constant, and angular 
eccentricity respectively. 

The main difference between the relativistic orbital representation and the non-rclativistic one is the appearance of 
three eccentricities e^, e^, and compared to one eccentricity in the Newtonian case. However, these eccentricities 
are related. The explicit expressions for the parameters n, k, a^, et, e^, e^, ft, gt, and in terms of the 2PN 
conserved energy and angular momentum per unit reduced mass are given by Eqs. (38) to (48) of @- Though the 
three eccentricities are related, there is the question of selecting a specific one to present the polarizations. We have 
chosen to present polarization waveforms in terms of et since it explicitly appears in the equation relating I to u, 
which we will numerically invert while computing its power spectrum. The question whether, as in pulsar timing, 
there is a particular combination of the three eccentricities 'a good eccentricity' in terms of which expressions take 
familiar Newtonian-like forms is interesting and open. 

Exactly as in the Newtonian case, at 2PN order too, one can split into two parts; a part A linearly increasing 
with time and a part W(Z) periodic in I but with a more complicated time variation: 

= A + W(O, (2.13a) 

A = 0o + (l + 4)^ (2-13b) 
k 

W(l) = (1 + — ) (v-u + et sinu) 

+ ^ |/0 sin2u -I- sinSv - ft sinv - gt{v - m)|. (2.13c) 

Note that com plica te d 2PN corrections to A and W(Z)also include k that represents the periastron advance. 

Using Eqs. (2.11), (2.12), and Eqs. (38) to (48) of [Q, it is straightforward to obtain the 2PN accurate expressions 
for the dynamical variables in terms of ^ = et and u, using the following relations, easily derivable from Eqs. 

(38) to (48) of g 
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-2E^ c^C^/^ jl + ^ 



-2Eh^ = (1 - e^) <{ 1 



15 — 1] 
120-48?/ 

^2/3 

^4/3 



^4/3 

'24~ 



15 — 15 ry — ry^ 



(2.14a) 



(17-7?7) 6(2 + 9 + 77 



- (360 - 144 T]) et Vl - + (225 - 277 + 29 77^) Ct'' 



24(l-e2)2 
- (210 - 190 r; + 30 77^) ej^ + 189 - 45 77 + 77 

3^2/3 ^4/3 



(1-e?) 4(l-e?)2 
eJl + e'/' (4-77) + 



(51 - 26 77)e2 + (78 - 28 77);>, 

^4/3 



(1152- 656 7; + 41 77^)64" 



^ 96(l-e2)3/2 

+ 1968 - 1088 7/ - 4 77^^ y/l - e2 + (720 - 288 77) (1 - e?) 

^2/3 ^4/3 



et 1 



2 24 (1-6^3/2 



(288 - 242 77 + 2l7;2)et2 



390-30877 + 217;^ i/l-e2 + (180- 72?/) {l-ef] 



(2.14b) 
(2.14c) 

(2.14d) 

(2.14e) 

_ dr 



Using Eqs. ( 2.14 ) in Eqs. ( 2.11 ) and ( 2.12 ), we obtain after some lengthy algebra, expressions for r, f — ^ 



1^ f , = A + W(/) and = f 1^ f , in terms of ^ = e*, .. given by: 



r = (^)V3(i_e, cosT.)^ 

71^ 

^4/3 



^2/3 



(6 — 7 77) et cos 7i + 18 — 2 77 



6 (1 — Ct cosw) 

(72 - 231 77 + 35 77^) (1 - 6(2) et cos w 



72 v/(l - e^)3 (1 - et cosu) 
-(72 + 75 7; + 8 7;2)et2 - 234 + 273 77 + 8 77^ ) ^/l - e^ 

-36 (1 - e^) (5 - 2 77) (2 + et cos u) 
I = u — et sin it 



^4/3 



1 



8 ^1 - e? (1 - et cosm) 
+12 (5 - 277) (m - v) (1 - et cosw)! , 
3^2/3 

W = f — 77 + et sin u + — 2V \ v — u + et sm u 

(1 ~ 6* 

^4/3 1 



et sin uJl-e^rj (4 + 77) 



(2.15a) 
(2.15b) 

(2.15c) 



4 W 1 — e( (1 — et cosw) 



(102 



32(l-e2)5/2 (i_et cosw)3 
-52 77)et2 - 156 + 56 77I et cos w + 77 (4 + 77) et^ + (l02 - 60 7; - 2 if) et^ 

+156 - 52 77 + 77^^ + (1 - ef) ^((3 e^ + 12) 77 - 8) (et cos 77)^ 

+ ((8 - 6 7;) et^ + 8 - 24 77) (e* cos it) - 12 77 et^ - (8 - 27 7?) et^ 
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et sin u + (1 — gf cos u)'^ 



;(l-e2) (5-2r;)-8 (51-26r;)et 



-78 - 28r; 



u + 8 {I — et cosu) 



(51-26r?) et^ 



78 - 28r; ./l-e2-(30-12r,)(l-2e^ + e^) 



(1 — et cosu) 

^4/3 1 



^2/3 

et sin-u<{ 1 H ^ (6-777) 



(72 - 231 7? + 35 77^) (et cos u)^ 



72 (1 — et cosw)'^ 
(216 - 693 77 + 105 ry2) (et cos u)^ + (324 + 513 77 - 96 t]'^) et cos 1 

(36 + 9 7;) 77 et^ - 468 - 15 7; + 35 77^ 



36 



(1 — et cos u) (4 — et cos u) (5 — 2 77) 
^2/3 



(1 — et cos 
^(4-77) et2 + 3 



1 



(1 — et) (1 — et cosu) 
1 

12 (1 — et cosu)^ 



(1 — 77) et cosu 
1 



(l-e?)3/2 



18 (1 — et cos-u) 



(et cosu - 2et2 + l) (5 - 2r;)| + L-^i ( -(9 - I977 - 14r72)et' 

J (1-ef) IV 



-36 + 2 77 - 8 77^ j (et cosu)^+ (^-(48 - I477 + 1777^)et 
+ (69 - 79 77 + 4772)et^ + 114 + 2?7 - 5 772^ (et cosu)^ 
+ ^-(6 - 32 77 - 772)et*+ (93 - 19 77 + 16 77^) e? 



222 + 50 77 + 7;^ (et cosm) - 6 77 (1 - 2 77) e^ 



(54 - 28 77 - 20 772)et*- (153 - 61 77 - 2 77^) et^ + 144 - 48 77| |. 



(2.15d) 



(2.15e) 



(2.15f) 



Note that the above 2PN accurate expressions for A and W(Z) in terms of ^, et, ..etc are expUcitly needed to explore 
the spectral decomposition of the polarization wavefo rms. They are not meant to be explicated in the 2PN expressions 
for h+ and h^, in terms of A and W(;) in Eq. ( 2.24 ) and ( 2.25|). 

The 2PN corrections to h+ and h^, in a form similar to Eqs. (2.10), are obtained using Eqs. (2.1). However, we need 
the 2PN corrections to hfj^ in j4_DM coordinates, as the parametric expressions for r, r, cj) and the = A + W(Z) split, 
give n by Eqs. (2.15), are in the ADM gauge. However, the 2PN corrections to hfj'", given by Eqs (5.3) and (5.4) of 
|p9| , are available only in harmonic (De-Donder) coordinates. Using, in a straightforward manner, the transformation 
equations of Damour and Schafer Ew to relate the dynamical variables in the harmonic and the ADM gauge, we 
obtain the 2PN accurate instantaneous contributions to hj^^ in the ADM gauge. For completeness, we list below the 
relevant transformation equations relating the harmonic (De-Donder) variables to the corresponding ADM ones. 



Gm 



7; + 2 (1 + 1277) 



Gm 



Vd = VA 



Gm 
Gmf 



(2.16a) 
(2.16b) 



7v^ + 38 



Gm 



3r' 



Gm 



8 



Gr. 



34- 



.Gm 



„ Gm] 



rD = f A 



Gm 
8c 



2 (1 + Urj) 19 77f 



(2.16c) 
(2.16d) 



The subscripts 'D' and 'A' denote quantities in the De-Donder ( harmonic ) and in the ADM coordinates respectively. 
Note that in all the above equations the differences between the two gauges are of 2PN order. As there is no difference 
between the harmonic and the ADM coordinates to IPN accuracy, no suffix is used in Eqs. (2.16) for the 2PN terms. 



Using Eqs. (2. 16), the 2PN corrections to hf in ADM coordinates can easily be obtained from Eqs. (5.3) and (5.4) of 



|p9| . For economy of presentation, we write (/i^ )a in the following manner, (/if, )a = ih^j)o+ ^Corrections'' , where 
(/i^-^)a represent the metric perturbations in the ADM coordinates. {hf^)o is a short hand notation for expressions on 
the r.h.s of Eqs. (5.3) and (5.4) of where N, n, v, w^, f, r are the ADM variables Na, ua, va, va, respectively. 
The ^Gorrections^ represent the differences at the 2PN order, that arise due to the change of the coordinate system, 
given by Eqs. ( ^.16| ). As the two coordinates are different only at the 2PN order, the ^Gorrections^ come only from 
the leading Newtonian terms in Eqs. (5.3) and (5.4) of ||39|. 



TT\ 

ij )o 



G Gm 



c^R2c^ TA 



(5 - 55 fi) 77 + 2 (1 + 12 77) 



Gm 



G m , rprp 



("-ij)A 



(14wi-6r-i)?7 + 8(l + 5 77) 
(10 4-18r-i)?7 -(4 + 68 7;) 



G m 



TA 

G m 

TA 



fA{n(iVj))A 



TT 



(2.17) 



To check the algebraic correctness of the above transformation, we compute the far-zone energy flux directly in the 
ADM coordinates using 



\~dt 



327rG 



(2.18) 



After a careful use of the transformation equations, the expression for {d£ / dt)A calculated above, matches with the 
expression for the far-zone energy ffux, Eq. (4.7a) of |p9] obtained earlier. This provides a useful check on the 
transformation from {hfj^)j^ to {hj^)j^. 

We now have all the inputs required to compute the 2PN corrections to h+ and in terms of a elegant and 
convenient parametrization using Eqs. (2.1). As mentioned in [0,R9|, there is no need to apply the TT projection to 



(hf ) given by Eq. (2.17) before contracting with p and q, as required by Eqs. (2.1). Thus, wc schematically write. 



lTT 



(2.19) 



The polarization states /i+ and h^, for Eqs. (2.19) are given by 



— ((P.v) 



1 



aw (P-v) (q.v) 



(q.n) 



(q.v)^)+^((p.n) 



— ( (p.n)(p.v) 



(q.n)(q.v) 



an 



. (P-n) (q.n) + — ( (p.n) (q.v) + (p.v) (q.n; 



(2.20a) 



(2.20b) 



Before proceeding to a lengthy but straightforward computation of the 'instantaneous' 2PN accurate polarizations 
h+ and hx, we anticipate the structure of the final result by schematically examining the functional forms in the 
intermediate steps of the above calculation. The polarizations in terms of has the form 



(•••) + (---)cos20+(---)sin20 



(•••) + (•••) cos 2(/) 



+ 



(•••) + (•■■) cos 20 -h(---) sin 20 
[(•••) cos 2(j>+{-- •) sin 20] + a„ 



(• • •)sin20 



I cos 20 -I- (• • •) sin 20 



(2.21a) 
(2.21b) 



9 



In the above and what follows (• • •) denotes a dependence on variables et,n, mi,TO2,i and u. The structure of the 
PN expansion of the coefficients above is the following: 



1 H — (• • •) cos(/) +(•■•) sir 



(• • •) cos 3(f>+ (• • •) sin 30 



(•••) + (•• •) cos 20+ (• • •) sin 20 
(•••) + (•••) cos 40 +(•••) sin 40 



+ 



(2.22) 



a„„ has a similar expansion. also has a similar expansion but the leading order term is of order - . Using this 
information about the functional dependence on and elementary trigonometry one can infer in detail the harmonics 
of a ppearing at each order and consequently the A and W(/) dependence on display in the equations (2.24) and 



( 2.25 ) below. The 'instantaneous' 2P N ac cu rate p ol arizat ions h + an d hx in terms of ^ = '^^^ ,m,ri,et and Sines and 



Cosines of X,W{l),u{l) and i, using ( |2.20D , (|1|), ( |2.15| ) and (|1|) are finally written as:' 



(2.23) 



where the curly brackets contain a post-Newtonian expansion. The expressions for various post-Newtonian terms in 
the 'plus' and 'cross' polarizations are shown below in a form emphasizing the harmonic content and the corresponding 
amplitude modulation. The various post-Newtonian corrections to the 'plus' polarization are given by: 



H 



(0) 



pO 

^C2C2 



cos2W + Pg2C2 sin2W 
sin2A + P°l , 



cos 2 A 



Pc2S2 cos2W + Pg,s, sin2W 



H 



(0.5) 



pofci cosW-KP^fbi sinW 



0.5 



cos A -I- 



P^ctsi cosW + pOfsi SinW 



sin A 



P^3^C3 cos3W + P^3^c3 sinSW 



cos 3 A - 



P^3^S3 cos3W + P°35s3 sinSW 



H 



(1) 



sin 3 A 

P^2C2 cos2W + P^2C2 sin2W 



-P^2S2 sin2W 



H 



(1.5) 



+P 



sin 2 A 

Pc4S4 cos4W-hPs4S4 sin 4 W 
Pi,fci cosW-HPi'i^ci SinW 
Pj,|c3 cos3W + P^35c3 sin3W 
sin 3 A 

P^|s5 cos5W + pi|s5 sinSW 
Pc2C2 cos2W + P^2C2 sin 2 W 



cos 2 A 

Pc4C4 cos4W-|-Pg4C4 sin 4 W 
sin4A-f P^ 



Pj,2S2 C0S2W 



cos A 



S3S3 sm3W 



hPi2S2 sin2W 



S2C2 

sin 2 A 



cos 4 A 



Pj,fsi cosW-l-P^fsi sinW 
cos 3 A 

P J;|c5 cos 5 W 4- Pl fc5 sin 5 W cos 5 A 
sin 5 A 
cos 2 A 



sin A 



Pj,|s3 cos3W 



Pl^S2 C0S2W 



Pc4C4 cos4W + Pg4C4 siu 4 W 



COS 4 A 



Pc4S4 COS 4 W Ps4g4 siu 4 W 



sin 4 A -t- 



P^6C6 C0S6W 



+Ps6C6 sin6W 



cos 6 A - 



P^ese cosGW-l-Piese sin6W 



sin 6 A + pH , 



(2.24a) 



(2.24b) 



(2.24c) 



(2.24d) 



(2.24e) 
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and for the 'cross' polarization by 



H 



(0) 



H 



(0.5) 



H 



(1.5) 



XO2C2 COS2W + X02C2 sin2W 



cos 2 A + 



Xc2S2 COS 2 W 



+Xs2S2 sin2W 



sin 2 A 



X^fci cosW + Xgfci sinW 



I yO.5 
+ ^C3( 

+X§35g3 sinSW 



cosSW + Xgg^cs sinSW 



cos A + 
cos 3 A 



X^fsi cosW + Xgfsi sinW 



sin A 



X^|g3 cos3W 



sin 3 A 



^(1)^ 



Xc2C2 cos 2 W + Xg2C2 ^in 2 W 



^C2S2 cos2W 



+Xg2S2 sin2W 



sin 2 A 

Xc4S4 cos 4 W + Xs4g4 sin 4 W 
Xj^fci cosW + X^fci sinW 



cos 2 A 

Xc4C4 COS 4 W + Xg4C4 sin 4 W 
sin4A + Xi 



sin A 



y1.5 
^C3C3 



COS A + 



cos3W + X^35c3 sinSW 



+X^3^g3 sinSW 



y1.5 

^C5C5 



sin 3 A 

Xi;fs5 cos5W + Xi55s5 sin5W 
X2 2C2 COS2W + X22C2 sin 2 W 
sin 2 A 



cos 4 A 



Xj^fsi cosW + X^fsi sinW 
cos 3 A 

cos 5 A 



Xcisa cos3W 



cos 5 W + y^f(^r.r. sin 5 W 
sin 5 A 
cos 2 A 



X2 2S2 cos2W + X22g2 sin 2 W 



Xc4C4 COS 4 W + Xs4C4 sin 4 W 



cos 4 A 



Xc4S4 cos4W 



+Xs4S4 sin4W 



sin 4 A 



X^gce cos6W + X2gc6 sinew 



cos 6 A 



X^ggg cosGW + X^ggg sin6W 



sin6A + X^ 



(2.25a) 



(2.25b) 



(2.25c) 



(2.25d) 



(2.25e) 



where P's and X's are functions of et, mi, m2, u{l) and i. The notation Pci^g;., for instance, denotes the coefficient of 
cosbW sine A at 'aPN' order and similar explanation holds for X^i^g^. too. The explicit expre ssion s for P's and X's, 
i.e , the coefficients of Sine and Cosine multiples of A and W(Z) appearing in Eqs. (2.24) and (2.25) are given by 



i^C2C2 - -i^S2S2 - (i^e^ eos^x)2 

pO — pO _ 9 1 

i^C2S2 - i^S2C2 - ^ cOS7.)2 



(1 + ^2) 



(et cos u) — (et cos u) — 2 + 2 



(1 + C^) \Jl-e1 et sinwl , 



52 



(1 — et cosw) 



etcos u , 



1 



po.5 _ po.5 _ -5* rni~m2 

^cici — ~^sisi — t; T\ V 

2 TO (1 — et cos u) 



(1 - 3C^)etsinu \ , 



(2.26a) 
(2.26b) 
(2.26c) 
(2.26d) 



pO.5 _ pO.5 



5* TOi — TO2 

4 



1 



TO (1 — et cos uy 



2 (cf cosm) + 5 



11 



-(6 (et cosu) - 1)C^ 
S mi — 1712 



2 m (1 — et cosu)^ 
-46*2 + 5 

S mi — 7712 

4 m (1 — et cosu)^ \ 



(1 + C'^)et sinw 



(et cos u) — 2 (et cos u) 



6 (et cosu) — 7(et cosw) 



-8e? + 9 



1 I / 

6 (1 — et cosu)^ 1^ 



— (9 — 5 77)(et cosu)^ + (18 — 10r?)(et cosu)^ 



+ ^(18 - 10 T])et^ + 20-12 T]^ (et cos u) - (33 + 11 T])et^ - 14 + 38 



+ 



-(3 + 13r?)(et cos«) + (6 + 26r?)(et cos«) + {6 + 26 r])et^ +33 



-51 r? (et cos u) - (48 - 34 r])et'^ + 6 - 22 r; 



C^ + {l-3r]) 



-6(etC0S'u) +12(etCosu) + (12et2 — 13)(et cosu) — Qet^ +4 



1 et sinu 



{A/r^[( 



6 (1 - et cosu)3 t y/l - e| 
-(39-7r?)et2 + 27 + 5r; 



-(30 + 20 r/)et2 + 18 + 32 r? 



(18 - 10 r7)et - 6 - 2 r/ (e* cos u) 



(6 + 26 r/)et2 + 6 - 38 r? (et cos u) 



C^ + 3^/l-et^l-3r1) 



-4 (et cosm) 



+5 



24 (1 - et 



^|(l-3r?) (1 + C^ 

cosu)* [ 



—6 (et cosu)* + 18 (et cosu)^ 



+(48 et^ - 61)(et cosu)^ + (65 - 69 et2)(et cosu) - 48 et^ + 117 et^ 



-64 



1 



4 {1 — et cos m)* 
+9 (et cosm) + 8et2 - 13 



(l + C^) (l-3?7)^l-e2 et si 



sinu 



—4 (et cosm) 



1 



24 (1 - etCosu)3 
- (57 + 5 T]) (et cos u) + (87 + 3 r])et^ 



(30 — 2r]){et cosu)^ — (60 — 4ry)(et cosu)^ 



+ 



(18- 54 77) (et cosu)^ - (36 



-108r7)(et cosm)^ + (3 - 9r7)(et cosu) + (15 - 45r?)et2 

S mi — m2 et sinu 



48 (et cosu)^ — 144 (et cosu)^ 



48 m (1 - et cosu)^ 
+(33 + 22 77)(e( cosm) - (336 + 12 7?)e(2 + 399 - 10 7? 



12 



+ -(108 + 72 7?)(et cosu)^ + (324 + 216r?)(et cosu) 
-(12 + 240 77)(et cosu) - (144 - 36 77)6*^ - 60 



+60 r? 



C^ + 5{l-et cosu)(l-2r7) 



pi. 5 _ pi. 5 _ 

^cisi — -t^sici — 



-24 (et cosm) + 5 

S mi — 7712 



96 



m {I - et 



12 {et cos u) 



(48 + 48 ry - 96 ej^ ) (e* cos w)" 



- (24 + 48 7? - 72 et^r?) (e* cos u) + ( (396 + 104 r])et^ - 204 - 296 ) (e* cos u) 

1 



+ 



(216 



+(123 - 198r?)e/ - (546 - 22Qrj)et^ + 303 + 98r? 
+144 7?)e(^ - 72 - 288 r/^ (e* cos m)^ - ^(756 + 240 j])et'^ - 444 - 552 77 
(et cosm)^ - ^(144 + 96r?)et^ - 336 + 96 r?^ (e* cosu) + (720 + 1447?)et 

120 (et cosuf 



- (756 + 96 r;)et^ + 12 -24 77 



l-e? (1-2,7) 



pi. 5 



C3C3 



pi. 5 



-276 (et cosu) + 52 (et cosm) + 105 e; - 1 
S mi — m2 et sin u 



S3S3 



96 



m (1 — et cosu)^ 



(108- 24r7)(et cosw)-" - (324 - 72 r;)(et cosu) 



(432 - 96r7)et^ - 143 - 2747? (e* cosu) + (1176 - 72r7)et^ - 671 



-346 77 



+ 



(72 + 48 77)(et cosu)^- (216 + 144r7)(et cosw)^ 



(288 



pi. 5 _ pi. 5 _ 
J^^C3S3 — J^^S3C3 — 



+192r7)et^ + 88- 7367? (et cosm) + (1152 - 24r7)et^ - 632 - 4247? 



+5 (1 - 277) 

+72 et^ - 77 

S mi — m2 
64 



12 (et cos u)^ - 36 (et cos uf + (77 - 48 e*^) (e* cos u) 



1 



m (1 - et cos 7^)4 [ y^i _ e2 



1 



(2.26m) 



(2.26n) 



(2.26o) 



^(216 - 4877)et^ - 120 -48 77^ (et costz)^ 



- ^(644 - 120 7/)et^ - 436 - 88 77^ (e* cos uf - ^(288 - 64 77)et'^ - (76 

+232 77)et^ - 340 + 4247?^ (e* cosm) + (719 + 277)et* - (510 + 43677)et^ 
1 



-225 + 450 77 



+ 



^(144 + 96 77)et^ - 48 - 192 77)(et cosTi)^ 



- ^(504 + 160 7;)et^ - 296 - 368 7/ j (e* cos uf - (^(192 + 128 77)64^ + (96 
-576 ??)et^ - 416 + 576 77 j (et cos u) + (928 - 416 77)et^ - (1016 



-576 77)6(2 + 72 - 144 7; 



^l-e| (1-277) 



120 (et cosu) 



13 




14 



+5 -(1800- 2856 77- 7128 ?7^)(et cosuf + (7200 - 11424 - 28512 r?^ 
){et cosu)*+ ((3600- 5712 77- 14256 r?2)et2- 4746 - 7342 7? + 61614 



(et cos uy + (^-(18684- 43500 77- 34452 77^)et - 6 + 23358 77 
-67518 77^^ {et cos uf + ((12909 - 32529 t] - 16815 ri'^)et^ + 8109 
-30609 77 + 19281 rjA (e* cos u) + (9045 - 24345 77 - 7335 r]^)et'^ 



(j4 



-(15915 - 43431 77 - 11289 77^)6*^ + 288 - 1184 77 + 672 77^ 



-360 (et cos u) + 1440 (e* cos u) + 720 e*' 



+ 15(1-577 + 5772) 



-2418Vet cosm)^ - (2124 6*^ - 2178) (e* cosu)^ + (1553 6*^ 



P2 "02 
C2S2 — -^^8202 



-527) (et cos u) + 345 et'^ - 839 e^^ + 32 
1 Cf sinii f 1 



(504 - 1048 77-40 77^)et* + (144 



192 (l-e^ cost/,)« I ^(1 - c^ f 
-592 77 + 48 77^)64^ - 1080 + 216 77 + 104 77^^ (e* cos uf + ((1584 
-2800 ?7 - 448 77^)64'' - (1536 - 1568 77 - 416 77^)64^ + 4272 - 784 77 
-544 77^^ (et cosuf+ ((1956 + 908 77- 92 77^)6*^ - (4584- 1736 77 

-760 r]^)et^ - 1692 - 628 77 - 92 tj^^ (et cos u) + (1359 - 6699 tj 
+171 77^)64^ - (7113 - 20941 77+ 13 77^)64^ + (10053 - 22809 77 



711 T]'^)et^ - 2859 + 7895 77 + 361 rj' 



2\„ 2 



(408 



-1464 77 - 776 r/^)ef - (1296 - 2896 77 - 2832 7/^)6*^ - 552 - 760 77 
-1864772](et cosuf + (-(1536 - 582477 - 200077^)64'' + (3648 



-11456 77 - 7648 77^)64^ + 2208 + 3616 77 + 5072 77^ (et cosm) 



+ (1692 + 2740 ?7 - 3012 77^)64^ - (3096 + 5768 77 - 9480 77^)64^ 



-2916 + 504477 - 5892 77^ J (e4 cosm) - (1167 - 98777 - 160577^)64^ 

+(2937 - 10061 77 - 3027 77^)64* - (2757 - 17289 77 - 151 ri'^)et'^ 

1 



+2427 -8887 7; + +1079 7;' 



(600 - 952 7/- 2376 77^)6t' 



-216 - 584 77 + 3528 77'' j (et cos u)"" + (2688 - 5440 77 - 7296 77^)64^ 
+1440- 44877- 11040 77^^ (64 cosu)^ + ((1684 - 25967/ - 6940 77^)64^ 



(2.26s) 



15 




16 



- (680 - 1936 r? - 632 r?^)et ^ + 224 - 336 r? - 1232 rj^ (e* cos u) 



-(879 - 2665 r/ + 93 r/^)e(^ + (1448 - 4084 r? - 788 ??^)et - 509 + 1179 r? 
+1061 ??2 + ( -(120 - 344 r/- 72 7/2)6(2 + 72 - 152 7; 

-216 7/2^ (et cosm)^+ ^(486- 14307?- 162 7/2)6*2 -330 + 8067? 
+630 7?2 ] {et cos uf + ( (240 - 688 7? - 144 T]^)et'^ - (540 - 1236 rj 



-1332 772)6(2 + 132+ 124??- 1692 7?2 \ {et cosu) - (936- 2950 7? 



+378 772)64* + (1566 - 4674 77 - 198 772)642 - 570 + 1484 77 + 756 772 



+v/l-et2(l- 5 77 + 5 7?2) 



72 (et cosu)^ — 240 (et cosu)2 



-(1446(2 - 364)(6t cosu) + 2496(2 - 301 



(1 + 6*2) (1 - 5 77 + 5 7?2) \ 120 (et cosm)^ 



1920 ' ' ' ' ' ' (1 - et cosu)6 

-600 (et cosu)'"^- (2160 et2 - 3206) (et cos ?i)'' + (7860 ef2 
-8444)(6( cosm)^ + (57606('' - 19135 62 + 13051)(6( cosm)2 - (114756(* 

-23240 6(2 + 11269)(6t cosu) - 3840 e*'' + 17235 6(* - 21325 6^ + 7776| , 
^'^^ 5^1 + C'2) (1 - 5 77 + 5 7?2) (72 (et cos uf 



192 ' ' ' ' {1-et cosm)6 

-312 (6( cosu)^ - (384 6(2 - 844) (6( cosu)2 + (1053 6(2 - 1325) (6( cosu) 



+384 6(* - 14376(2 + 1105 



W2 (l-e.'cos„). { [(120 - 120 , - 8 ,^)(., c„s„)'' + (3360 
-2608 77 + 416 772)(6( cosw,)'' - (11058 - 7718 77 + 1086 7/2) (6( cosw) 

2 



+ 



^(168 + 1240 77 + 152 7/2)6(2 + 13710 - 8746 7/ + 738 7/2^ {et cos u) 
+ ^-(3441 - 139 7? - 277 r]^)et^ - 5181 + 2751 7? - 423 rj^^ {et cos u) 
+(135 - 47777 - 363772)et^ + (3003 - 425 7/ + 2977/2)6(2 - 816 + 528 77 
+2 (180 - 516 7/- 108 7/2) (e( cos m)^ - (720 - 20647/ - 432 7/2) (e( cosu)* 
+(678 - 192877 - 504772)(6t cosu)^ + ^(150 - 430 7/ - 90 7/2)6(2 + 336 
-1010 77 + 90 7/2^ {et cosm)2+ ^-(930 - 2705 77 - 315 7/2)6(2 
-96 + 283 7/ + 9 7/2 ] {et cos u) + (378 - 11077/ - 81 7/2)6(* + (24 - 61 7/ 



-63 7?2)6(^ 



C2 + (1-577 + 5772) 



120 (6( cosTi) - 480 (6( cosu) 



17 



+566 (et cos u)^ + 84 - 102 (e* cos m) - 343 - 1 (e* cos u) 



12 (l-etCosu)2\(l-e2) [V 



+105 6^ + 49 6*2 

(1 - 6?) (240 - 193?7 + 24?72)j (6* cosu) + 51 - 33r? 
-(60- 24r?) 6t COSM + 150- 60?? 



+ 



6* COSM 



^C2C2 



^0 



■S2S2 



C2S2 



S2C2 



yO.5 



ClCl 



yO.5 
^ClSl 



yO.5 
- A 



SlSl 



yO.5 
^SlCl 



yO.5 



C3C3 



^C3S3 



^C2C2 — 



Xc2S2 



yO.5 
- A 



S3S3 



yO.5 
^S3C3 



"XS2S2 



^S2C2 



C4C4 



y1 - 

^C4S4 — 



■S4S4 



^S4C4 



-4C 

2C 



(1 — 6t COSU) 
1 



(1 — et cosu)2 
(75 mi — m2 a/I 



2 m (1 — 6* cosu)^ 

, mi — 1712 1 



2 6t sm u , 
{et cosu)2 — {et cosu) — 2 (et^ ~ "'"^} 
2 (6t cosw) — 3 > , 



-CS- 





m 


cs 


mi - 


- m2 




m 


-8 6*2 


+ 9 


}• 


-CS 


mi - 


- m2 


m 


-W 


+ 5 




c 


6* sin It 



6* Sin u , 



6 {et cos u) — 7 (6t cos w) 



3 < e* sm u 



{et cosw) — 2 (e* cosu) 



3 (l-6t cosw)3 lyr^ 
+ (33 + 11 r?)et2- 21 - 23 7/ 



-(12 + 877)642 + 2077 (6t COSU 



C 



6 (1 — 64 cosu)" 



+ 3^/1-6^ (1-377) (2(6t cosw)-3)C2 

,3 , tnA , 1^ \l- \2 



-(12 + 8 r?)(6t cosu)'' + (24+ 167?)(6t COSM) 



+ (24 + 16 77)64 + 53 - 63 77 (e* cos u) - (69 + 13 77)6* 



-20 + 5277 +(1-377) 

-(13 - 12 6^) {et cos w) - 21 6^ + 16 
6t sinu 



6 {et cosu)^ + 12 {et cosu)^ 
C' 



2 (l-etCosu)4 
-8 6*2 + 13|, 
CS'2 1 



1 — 62 (1 — 3 7;) < 4 (et cos u) — 9 {et cos w) 



(1 — 3r7)< —6 {et cos it) 



12 (l-etC0S'u)4 
+18 {et cosu)^ + (48 6*2 — 61) (6t cosu)^ — (69 6(2 _ 55^ costx) 

-48 6t^ + 1176*2 -64|, 



18 



y1-5 
"^SlSl 



X, 



1.5 
SlCl 



y1-5 
- A 



S3S3 



X 



- A 



S5S5 



y1.5 
^S5C5 



Cf sin u 



CS^ 

2 (1 — et cosu)^ 
C S mi — 1712 1 
Is rrT 



l-e2(l-377), 



-(96 6*2 -48- 48 r?)(et cosm)^ 



(1-et cosu)4 {^i- e| 
+ ( (432 - 24 r?)et2 - 264 - 144 r]] {et cos m)^ - ( (84 + 88 r])et^ 



+108- 280 7? j(et cosm) - (477 - 138 r/)et^ + (702 - 164??) ct^ 

-153 - 467? + \/l - 6*2(1 -2r]) 24 (e* cosw)^ - 84 (e* cosw)^ 

48 (et cosm)^ - 144 (et cosm) + (33 
+ (1-2??) 



+68 (et cos u) — 3 et 2 — 5 
C S mi — 1712 Ct sin u 



24 m (1 - et cosu)^ 
+22 r/)(et cosm) + (216 + 36 r/)et2 - 153- 58 r? 



12 (et cosu)^ — 36 (et cosu)^ + 29 (et cosm) + 24 et^ — 29 



C S TOi — m2 
"32 m 



1^(168 + 48 r?)et2 - 72 



(1 - et cosm)4 1^ y/l 
-144??^ (et cosu)^ + ^-(540 + 88r?)et2 + 332 + 296 r?^ (e* cosu) 

+ (- (224 + 64 77)et'' - (60 - 504 Tj)et'^ + 412 - 568 t]) {et cos u) 

+(725 - 10 r/)et* - (570 + 316 r])et^ - 171 + 342 r? 

72 (et cosu)^ + 172 (et cosu) 

C2 



1.5 _ 
S3C3 — 



'1-e? (1-277) 

-140j (et cosu) - 191 et^ + 135 

C S mi — m2 et sin u 
48 



96 ef 



(84 + 24 r/)(et cosit)^ - ^252 



m (1 — et cosu)^ 
+72 7?^(et cosuf - ^(336 + 96r?)ef2 + ll-582r?^(et cosu) 



+ (1080 + 120 r])et^ - 565 - 558 r? 



+ 3(l-2?7) 



12 (et cosu)^ 



-36 (et cos u) - (48 et^ - 77) (et cos m) + 88 ej^ - 93 



c- 



CS^ 



mi — m2 



{l-2r,) 



96 m (1 — et cosu)^ 

120 (et cosu)"- 396(et cosw)^ - (480 et'^ - 808) (e* cosu)^ + ( 825 et^ 



-773 J (et cosu) + 384et* - 1113 e*^ + 625 
C mi — m2 et sin u 



(1 - 2r7)<^ 12 (et cosm) 
48 m (1 - et cosm)^ ^ V t 7 

-48 (et cosu)^ - (144 et^ - 209) (et cosu)^ + (360 et^ - 394) (e* cosw) 
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Xc2C2 — ~^S2S2 — 2 C 



+192et^-600et^+413 
et sinu 



(1 — et cosu)^ 

C ei sin u 



[ 5-2r? 



(2 6*2 + l)(et cos m) +8 6*2 - 5 
(24 - 568 7/ - 8 r/^)et^ 



96 (1-et cosu)5 l(l-e|)i 
+(720 - 1872 r? - 720 7?2)et2 + 744 + 632 + 520 r?^ ) (et cos u)^ 



+ 1^(336 - 3056 Tj + 352 7]"^)^^ - (1728 - 7840 77 

- 1504 772)64^ - 2928 - 2768 r? - 1280 ??^^ (et cos u)^ + ^- (2652 
-988 7] - 852 ?72)et^ + (5400 - 3224 rj - 4008 r]'^)et'^ 
+1572 + 220 ?7 + 2580 r?^ j (e,, cos u) - (885 - 5145 r) + 297 ?72)et^ 
+ (4995 - 13935 ?? - 321 r/2)et'^ - (7047 - 12691 rj 



+ 



-(216 



-2333 ry2)et^ + 1497 - 3229 ry - 1523 77^ 



-568 7] - 264 r7^)et^ + 120 - 184 - 552 rj^ {et cos uf + (936 



-2552 ?7- 840 r?2)et2- 600 + 1208 ?7+ 1848 ?7M(et cosu)^ + -(868 



-2220 r] - 1220 r]^)et^ + 484 - 684 r] - 2372 77^ j (ej cos m) - (597 
-1737 77 - 303 77^)6*'* + (1342 - 3710 77 - 1250 77^)6*2 - 601 



+1397 77+137977^ 



C2 + 0T;2(1 - 5 77 + 57?2) 



120 (et cosu) 



432 (et cosu)^ + 484 (et cosu) + 75 et^ - 247 

c 1 r 1 



X^2S2 - Xi2C2 - ^ (l_et C0S«)5 I (1^) 



(1-e;) (360- 8520 77 

120 77^) j(et cosw)^- ^(1 - e?) (1440 - 34080 r; - 480 7?^) 

(et cosu)^ + ^(720 - 17040 - 2mrf)et^ 
-(43158 - 7891077 - 46290 772)et^ - 31002 - 14350 ?7 

-46050 77^^(64 cosu)'V ^-(16020 - 146340 77 + 21540 772)et'' 
+ (119994 - 295890 rj - 79710 Tf)et^ + 116346 + 6990 77 
+101250 7/2^ (et cosu)^ + ^(110835 - 10935 77 - 18105 77^)64^ 
- (276384 + 12480 77 - 145680 772)et^ - 54771 + 165975 77 
-127575 77^^ (et cosu) + (26955 - 16177577 + 45135 7/2)et*' 
-(176400 - 366960 77 + 95520 77^)64'' + (279333 
-230305 77 + 23505 772)6*^ - 56448 - 22400 77 + 26880 77^ 



(2.27q) 



(2.27r) 
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- ( 1 - e? ) (28800 - 1 1520 ?7) (et cos u) 



+ (^{1 - e^) (158400- 63360 7?)^ (et cosw)^ - |^(57600 - 23040 7y)et* 

-(244800 - 97920 77)et^ + 187200 - 74880 r?^ (e* cosit)^ 

+ ^(345600 - 138240 ?7)et^ - (417600 - 167040 77)64^ 

+72000 - 28800 r/^ (e* cosw)^ + ^-(518400 - 207360 ?7)et'' + (590400 

-236160 r?)et^ - 72000 + 28800 77 ] (e* cos u) + (230400 - 92160 r])et^ 



-^3240- 8520 7/ 



-3960 ?7^j {et cosuy + (12960 - 340807? - 15840 77^) (e* cosu)* 
+ ^(6480 - 17040 7? - 7920 r]^)et^ - 12582 + 24070 77 + 43770 77^ 
{et cos uf + (- (35820 - 101340 77- 21060 772)64^ +6846 



-(288000- 115200 7/)et2 +57600- 23040 77 



+4530 77 - 78930 77^^ (e* cos uf + ^(28455 - 74415 77 - 33825 772)6*^ 

+12159 - 64695 77 + 85575 77^^ (ct cos u) + (20655 - 68535 77 

+21735 77^)64'' - (40425 - 127185 r/ + 22785 r]^)et'^ + 4512 - 6880 r/ 



-16800 77^ 



C^ + 15 (1-577 + 577^) 



-120 (et cosu) + 480 (e* costx) 



+(240et2- 806)(et cosu f- (900 6(2 -918) (et cosw) 



X 



C4C4 



— — 

— ~^S4S4 — 



+(955 et^ - 613)(et cosu) - 45 e** - 205 et^ + 96 

CS^ etsinii f 1 rA^„„ 

-;^7-. — , (120-34477 

24 (1-et cosu)Hyi^ LV 

-72772)64^-72 + 15277+216772^(6* cosuf + ^-(480- 1400 77 

- 192 77^)6*2 + 324 - 776 77 - 660 772^ {et cos uf + (- (240 - 688 77 
-144 772)64^ + (518- 112677- 1442 772)6*2 - 110-23477 
+1802772^ (64 cosm) + (831 - 242577 - l^lri^)et^ - (1340 

-3544 77 - 1328 772)642 + 449 - 879 7? - 1361 7?2 + \j\-e\ (1 - 5 77 



+ 5 7/2) 



^C4S4 — ^S4C4 



-2016*2 + 241 
CS"^ 1 



48 (et cosu)^ + 162 (ej cos 'ti)^ + (gOej^ - 250) (64 cosu) 
C2~ 



720 (1 - 64 cosu)5 



(900 - 2580 77 - 540 7/2) (et cosu)^ 



(2.27s) 



(2.27t) 
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^C6C6 — ~^S6S6 



v2 _ y2 
^C6S6 — ^S6C6 



(3600- 10320?/- 2160 r;2)(et cosu) 



(7200- 20640 ry 



-4320 ?7^)et^ +4830- 7240 ?7- 24120 772^(et cos uf + (^(27990 



-82350 rj - 7290 T]^)et^ - 12120 + 21390 ?7 + 47610 tj^ j {et cos u) 

+ ^(7200 - 20640 77 - 4320 T]^)et^ - (8430 - 12105 77 - 42675 T]^)et'^ 

-13500 + 61875 77 - 65535 77^^ (e* cos u) - (24930 - 74325 77 
+3465 ?7^)et'* + (23100 - 57765 77 - 24135 ri^)et^ + 5760 - 30080 77 



+32640 77" 



+ 3(1-577 + 577^) 



120 (et cosm)^ - 480 (e* cosu)" 



-(1526- 960 6*2) (et ^qsu)^ + (2700 e? - 2718) (et cosw)^ + (960 et" 



-3235 et^ + 1933)(ef cosm) - 2115 et"* + 3805 et^ - 1536 



9 X / o et smw , ,4 

^ (1 + 5 77^ - 5 77) ^1 - (^^^^^ I -72 (et cos .)" 

+312 (et cos u)^ + (384 - 844) (e* cos uf - (1053 et^ 
-1325)(et cosu) - 384 ef-* + 1437 et^ - 1105| , 

^ (1 - 577 + 5 77^) (120 (et cos 7.)^ 

960 (1 — et cosit)° 1^ 

-600 (et cosu)^ - (2160et2 - 3206)(ef cosuf + (7860et2 

-8444)(et cosu)^ + (5760 et"* - 19135 6(2 + 13051)(6f cosu)^ 

-(11475 et"* - 23240 et^ + 11269)(et cosu) - 3840 et^ + 17235 et" 



-21325 et2 + 7776|, 



CS^ 



et smu 



24 (1 - et cosu)5 [ 



-(36-11677+12772)6*2 



+24 - 68 7/ - 24 77^ (et cos + (174 - 538 77 + 10 7j^)e/ - 150 



+442 77 + 62 77^ (et cos ?/) - (153 - 459 77 - 21 77^)et'* + (168 - 496 77 



-40 772)et2 - 27 + 85 77 - 17772 



+ \/l-e?(l + 5 7/2 - 5 77) 



6 (et cos u) — 22 (et cos u) + 15 et^ + 1 



where C = cosi and S — sini. The relations between the coefficients Hke PcnCn — ^SnSn -^CnSn 



(2.27u) 



(2.27v) 



(2.27w) 



(2.27x) 



SnCn 



are 



a trivial trigonometric consequence of the <j) — X + W split in the expression for GW polarizations in terms of 1 

To compare with the earlier 2PN accurate gauge independent expressions for hx and for binaries in circular 
orbits, we proceed as follows. First, we set et = in Eqs. (2.23) and rewrite the resulting expressions for h-^- and 
/ix in terms of the 'gauge independent' orbital angular frequency uj for circular orbits. The 2PN accurate relation 
connecting the mean motion n to uj may be derived from Eqs. (39), (44) and (46) of j44| and it reads: 



e = T 1- 



3 T3 



(l-e2) 4(1 -e2) 



(51 - 2677) 



(1-e?) 



(69 - 547?) 



(2.28) 



where r — . Next, we use the following angular transformation relation cj) = X + W{1) ~ (^biww ^ f ), where 
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'?i'Biww the orbital phase variable appearing in The expressions for hx and h+ thus obtained agree with Eqs. 
(2), (3) and (4) of |] modulo the tail terms. 
All the computations to obtain Eqs. ( 2.23| ) are performed using MAPLE |^^. This completes the calculation of the 



2PN accurate GW polarizations for compact binaries moving on elliptic orbits, modulo the tail terms 0. Though in 
principle the required equations for the tails are available in , the explicit expressions for the tail contribution to 
hy and for eccentric binaries have not been obtained. As mentioned earlier, this should be computed and included 
to write down the complete 2PN polarizations. 

III. INFLUENCE OF THE ORBITAL PARAMETERS ON THE WAVEFORM 

In this section, we investigate the dominant effects of eccentricity, orbital inclination and other orbital elements 
on hxit) and h^(t). For this purpose, the one sided power spectral density of the Newtonian contributions to the 
polarization waveforms are computed, by taking the squared-modulus of their respective discrete Fourier transforms, 
sampled over an orbital period. The results thus describe the influence of orbital elements on the power spectrum 
of Newtonian waveforms when gravitational radiation-reaction is negligible and referred to here as a 'non-evolving' 
waveform. 

To relate earlier studies done at Newtonian order to the present one, we proceed in two stages. In the first instance, 
to compare with the results of p9| , ^ , the orbital motion is restricted to the leading Newtonian order, and the 
periastron advance is mimicked by the introduction of an arbitrary constant shift parameter k in the 4> variable. In 
the second case, the orbital motion is taken to be 2PN accurate. In this case, the periastron advance is fully included 
in the formalism and explicitly defined in terms of the binary's parameters like the masses and eccentricity. In both 
the cases mentioned above, only the leading Newtonian part of the GW polarizations is considered. 

Let us begin with the 'x' polarization. For the ease of presentation, , the Newtonian part of hy {t), is written 
compactly below as: 

Grnn (0)^ (3.1a) 

ffx°^ = |a2s(0 sin2A + A2c(0 cos2a|, (3.1b) 

where A2s(0 = Msiu{l)) = X°,2S2 cos2W(/) + X§2g2 sin2W(0, and A2c(0 = Mc{u{l)) = X^.^^^ cos2W(/) + 
Xg2Q2 sin2W(^). Note that A2s(0 ^^'^ ^2c{l) are real and periodic functions of I with period given by The 
spectral analysis of hx will be performed using H^x \ scaled x polarization waveform, because since we are 
dealing with non- evolving binaries, i^^)^/^ essentially remains a constant over a few orbital periods. Similar 

arguments hold for h-^. too. 

A. Newtonian orbital motion 

In this section, we restrict the dynamics of the binary to Newtonian order. This implies we are using Eqs. (|2.8|) 



for A and W(^) in the 4> — X + W(/) split. However, following |29|32], we introduce an arbitrary periastron advance 
parameter k into the definition of A so that \ — <j>o + {1 + k) I and W(^) = {v — u + e smu) (1 -|- k). Note that with 
these forms for A and W(^), the scaled GW polarization waveforms are entirely specified by e, k and ^o- 

As mentioned earlier, A2s(0 ^'^d A2c(0 periodic functions of I = n(t — tg), where n = 27r/r, fr being the 
frequency associated with the 'radial period' p2| i.e. the time of return to the periastron. Consequently, they can be 
expanded in Fourier series as follows 



oo 



A2c(0= E ^J-^'''' (3.2a) 

OO 

A2s(0= E ^J^'''- (3.2b) 



j = -oo 



A C or Fortran version of the above hx and fi+ expressions is available on request from gopuQwugrav. wustl.edu 
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Employing, Eqs. (U) and X = (t)o + {1 + k) I, in (|^), we get 



(3.3) 



where 



Eq. (3.3) may be re-written as 



J2C 



C, EE 



{Cj — i Sj) , 



u;+ = {j + 2p), 
ujj = {]-2p), 
p={\ + k). 



(3.4a) 

(3.4b) 

(3.4c) 
(3.4d) 
(3.4e) 



E 

J=0 



S'^j e" 



I < iuj . I 



(3.5) 



Recalling I = nit — to), with n = 2t: fr, the frequency content and the associated intensities may be read off from 
the above. From Eqs.(t3.5|), it follows that the Fourier spectrum of H^^\l), consists of lines at frequencies uj^ fr and 



fr with powers {\SjY' + \C-j\ ) and (|5- 



\Cj\ ) respectively. Using the reality of (/)o, A2S and A2C, though Sj 



and Cj are complex numbers, it is easy to show that 



w+ fr will be 2 l^jf . Similarly, \Cj 

Thus, the Fourier series for Newtonian part of ft,x effectively reduces to 



|C_j| implying that power in the line with frequency 



-J I % and power in the uj^ fr line is 2 



4")=x/2rE 



D-i e J 



Cl UJ ■ 



(3.6) 



The 'one sided power spectrum' for the Newtonian /ix may be written as 



H^'^^^lSoe 



00 r 

E 



{i+2 P) I 



+ C, e 



i\U-2p)\l 



where explicitly the sum is over positive frequencies. In the generic case, the lu. 



frequencies ( 1 + 2 p) /j. , (2 + 2p) fr,{3 + 2p) fr, with strengths ~ | •S'l p , 

LuJ' = [j — 2p) part of Eq. creates hues at frequencies (1 — 2p) fr, (2 



\S^?, 

~2p)fr,{3 



(3.7) 

(j + 2p) part gives hues at 
... respectively. Similarly, the 
- 2p) fr, ....etc with strengths 



\Ci\ 



.etc respectively. There will be also a line at frequency 2pfr with strength 



proportional to 

These observations are easy to understand. At Newtonian order, in the absence of the periastron precession, i.e. 
k — 0, there is only one time scale in the problem, given by the orbital period and the spectrum consists of lines at 
multiples of the orbital frequency. When periastron precession is introduced, fc ^ 0, a second slower time-scale enters 
the problem, which splits and shifts original spectral lines from their earlier positions, thereby lifting the degeneracy 
associated with the non-precessing orbit. 



A caveat is worth noting: The discussion after Eq. (B.7) is valid only if all the terms corresponding to frequencies 



{js + 2p) fr and {\jc — 2p\) fr are linearly independent, where jc and js are summation index j for Cj and Sj. This 
is in general true except when js + 2p = \jc — 2p\, which corresponds to values of k equal to 0, 0.25 and 0.5 [^Oj. For 
these values of k, power in a given spectral line will have contributions both from Cj and Sj for different values of 
j given by jc — — 4 = 4 fc. These special values of k are interesting in that they can provide useful checks on the 
numerical accuracy of the analytical procedure outlined above. This is because for values of A: = 0, 0.5, 0.25, the full 
time domain waveform [ and not just the parts A2s(0; A2c(0 ] -^x (0 is exactly periodic over 2tt, 2n and Att intervals 
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respectively. Consequently, one may alternatively compute the desired power spectrum by a direct Fourier transform 
of the full h'^\i), without going via Eq. ( |3.7| ), which exploits double periodicity of 4> in I and A. Similar arguments 
hold true for the + polarization. 

It is clear from Eq. (^^) that the strengths of the different Fourier components are determined by the coefficients 
Cj and Sj which are given in terms of Cj and Sj, the discrete Fourier transforms of A2c(0 ^'^d A2s(0 0- This has 
become possible since we have exploited the double-periodicity of the motion in angles I and A. Thus the calculation 
reduces to the numerical implementation of Sj and Cj which we turn to next. 



Though the power spectrum for the Newtonian part of /ix can be obtained using Eq. (3.7), its implementation is 
not straightforward due to following reasons. First, the Discrete Fourier Transforms Cj and Sj can be evaluated using 
standard Fast Fourier Transform routines as in Numerical Recipes [ pT| only after A2S and A2C are written as explicit 
functions of I. However, in our analysis they are explicit functions o f u a nd thus implicit functions of I via I = u—e sinit. 
Consequently, we must first compute u(l) and substitute it in Eq. (3.1) to proceed. Secondly, W(m) = {v — u + e sin it) 

will not be a smooth function of u if we numerically implement v — 2 tan^^ ^ (i^f ) tan(|^)|. We will also need to 

use a smooth functional relation connecting v and u to obtain a well behaved -ff^ (""(0 ^^'^ H^{u{l). 

Let us first consider the implementation of u(l). There are two independent ways to obtain u(l) from I = u — e sinw. 
The first method is widely used, for analytical treatments, in standard textbooks of celestial mechanics |52|. The idea 
here is to expand the eccentric anomaly u in terms of the mean anomaly /. At the Newtonian order, it is given by 



/ + 



s=l 



Js{s e) sin si , 



(3.8) 



where Js{s e) is the Bessel function of the first kind of order s with s > 1. 

Alternatively, we can numerically invert Eq.(2.5) connecting the mean and eccentric anomalies, using the Newton- 
Raphson method implemented by rtsafe routine of Numerical Recipes pT[ |, and obtain u{l). We compute u{l) using 
both methods to make sure that they give consistent results, for the parameter values we are dealing here. 

We now turn to the numerical implementation of W(u). In text books of celestial mechanics, the transcendental 
relation connecting true and eccentric anomalies is expressed as a series given by 



u = 2 



— smj u 



(3.9) 



e sinu, to circumvent artificial 



where /3 — ^ (l — \/l — e^) . We use above expansion of w — u in W(m) — v ~ 
discontinuities in W(ii) as a function of the eccentric anomaly u. 

Using the above inputs, we compute u(l) and W(u(l)) at a finite number of points by sampling I. Next, we use 
the realft routine of ]5l| ] to compute the discrete Fourier transforms Sj and Cj of the discretely sampled peri odic 
functions A2s(0 ^-^d A2c(0- We then compute the 'one sided power spectrum' for the Newtonian hx using Eq. (3.7) 
for various values of e, fc and i. We now have all the inputs to investigate the influence of orbital elements on the 
Newtonian part of the x polarization waveform. The results and discussions are postponed to the end of this section. 

The spectral analysis for /i+ is similar to that for hx and we only quote the main results without any further details. 



H^^^ = Po(0 +Pi(0 sin2A + P2(0 cos2A, 

where Pi = P^2S2 cos 2 W + Ps2S2 sin 2 W, P2 = Pc2C2 cos 2 W + Ps2C2 sin 2 W and Pq 
e cosu). The Fourier series for Eq. (3.10) is given by 



(3.10) 

(1 - C2)e cosu(l - 



H 



(0) 



j = -OD 



S+e'^^ ' + C+e'^^' +P°e''' 



(3.11) 



where Sj' and C^' are defined similar to Sj and Cj but with ^2S and ^420 replaced by Pi and P2. Similarly, P? is 
the discrete Fourier transform of Po(0- Using arguments similar to the ones used for the analysis, we relate Sj', 



CTj and STj, 



C and obtain the 'one sided power spectrum' for Newtonian h+ as: 



* Using the Fourier integral theorem, it is easy to show that A2c(0 may be written in terms of the Fourier transform of 
^2c{l), the discretized version of which allows us to express A2c(0 i^i terms its discrete Fourier transform. Similar arguments 
apply to A2s(0- 
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Ij^i L J J 

From Eq. ( ^.12[ ), it follows that for the '+' polarization there will be lines at frequencies 0, 2p/r, |1 — 
Mfr, fr_,{l + 2p)U,_ |2 - 2p\U, 2/,, (2 + 2p)U. ... with relative strengths ^ ^\P^\\ \S+\', \Ct?, 

IC'^P, |-P2pj I ■5*2^ Pi respectively. Note that there are lines unaffected by introduction of k. These arise 

from the non-A term in H^{1). The values of fc = 0,0.25,0.5 are special and require a treatment analogous to the 
corresponding one in the cross polarization case. 

Using the above inputs, we plot the time-domain waveforms and the associated normalized relative power 

spectrum ^ (h'-'-''' y Figs. @-(|6|). The combined influence of the orbital parameters like eccentricity e, periastron 

advance parameter k and orbital inclination i, on the time-domain waveform and the associated power spectrum of 
the Newtonian hy and using Newtonian accurate orbital motion is summarized below. 

• Eccentricity, e . 

The effects of e on iJ" arc explored in Fig. (||). In the limit of low values of eccentricity e, as expected, the 
dominant contribution to the power spectrum, comes from the second harmonic. However, as the eccentricity 
e increases, higher harmonics appear in the spectrum with comparable strengths. For a given value of the 
periastron advance k and inclination angle i, the position of the dominant harmonic changes as the value of e 
increases. The shape of the waveform also changes significantly as we increase e. For moderate and high values 
of e there is a stronger burst of radiation near values I — and 2 tt corresponding to the periastron passage, 
since near the periastron, the two masses are closest to each other and their relative velocity is a maximum. 
In the frequency domain this results in the broad peak containing many frequencies. The line feature in the 
frequency domain on the other hand corresponds to the average orbital motion of the binary. 

• The 'arbitrary' periastron advance parameter, fc. The observation made here are based on Fig. (^). A 

careful inspection of ^(/) with fc ^ indicates that in general they are not 2tt periodic. This is expected as 
fc is a measure of the angle of return to the periastron. As mentioned earlier, in the power spectrum, the main 
effect of including an arbitrary fc is a 'splitting' and subsequent 'shifting' of the position of each spectral line 
from its integer multiple value in units of fr, the radial frequency. The shift is appreciable for medium and high 
eccentricities and leads to a shift of the dominant harmonic in the spectrum. 

• Orbital inclination, i. A change in the orbital inclination changes only the magnitude of and its power 
spectrum, keeping the relative distribution of spectral lines the same. This is easy to see as the dependence 
of orbital inclination angle i is easily factored out in the expression for hx- However, the shape of TJ^ and its 
power spectrum is influenced by i as seen in Fig. (^). 

If the polarizations of the gravitational wave hy and are available, the orbital inclination can be inferred by 

computing the ratio of the total power measured in each polarization. For circular orbits, the result analytically 

< 2 ■ 

follows since jj^w^ is only a function of i and is given by -pqr^^V^- In the general eccentric case, to explore this. 
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we plot jpy-p &s a function of i for various values of the e and fc in Fig. (|^). The plots are identical for different 

values of e and fc but vary with the inclination angle providing support for the claim made in the beginning of 
this paragraph. 

Even though in general, an arbitrary periastron advance parameter fc at Newtonian order destroys the 27r-periodicity 
of 77°, + (0 1 we may choose fc values exactly equal to or 0.5 so that (1) is still 2 n periodic. These particular values 
of fc allow us to perform useful numerical checks on our analytical procedure. In this case, we can compute the power 



spectrum directly from Hy{l) by numerically implementing the Discrete Fourier Transform of Eqs. (3.1). We can 



also implement Eq. (3.7) to obtain the power spectrum, after adding contributions from various Cj and Sj to a given 
harmonic, which is now some integer multiple of the radial frequency, fr- The results are displayed in Fig. (^). For 
better comparison, in these flgures, we normalize relative to the power in the dominant harmonic rather than relative 
to the total power as in other figures. We next, choose fc = 0.25, so that now ( l), is now Air periodic. Again, 



comparison with the power spectrum computed directly from Eq. (3.1) and via Eq. ( |3.7[ ) is possible. The results for 



fc = 0, 0.25 and 0.5 via these two methods are compared and found to be identical up to numerical errors as seen in 



Fig. 



), providing important checks on our analysis and routines that compute the one sided power spectrum via 



Eq.(p77). We observe a similar behavior for H^. 
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B. The 2PN accurate orbital description 



The spectral analysis discussed in the previous section may be extended to 2PN accurate orbital motion with 
minor technical modifications. The expressions for W(^) and A in the cf) = X + W(Z) split are now given by Eqs. 

( ^I3| ). Mor cover, the orbital elements appearing in 4")and<' are now 2PN accurate. These changes will modify 
expressions for Sj, Cj in Eq. (3/7) for the 'x' waveform and the corresponding expressions for the '+' waveform. 

To implement the 2PN accurate spectral analysis, we note the following: First, at the 2PN level the simpler 
approach to ob tain the u{l) relation, connecting the mean and eccentric anomalies, is to numerically solve for u{l) 



from Eq. (2.11) because 2PN accurate analytic expression for u{l) similar to Eq. (3.8) is not available in the literature 



However, we may employ Eq.(3^) with in the place e to get v — u ai 2PN order. This is because in the generalized 
quasi-Keplerian representation, the relation connecting true anomaly v to eccentric anomaly u has the same structural 
form as for the Keplerian case. Secondly, there are post-Newtonian corrections to the relations connecting Cr and 
to et, TOi, 7712 and 77. In our analysis, only those values of Ct are considered which lead to e^, less than one. Finally, 
in this 2PN accurate orbital description, the periastron precession constant k is no longer arbitrary but uniquely 
determined by 7771, 7772, 77 and e* as given by in Eq. (2.14). 

We explore the effects of 2PN accurate orbital motion on the power spectrum for 'x' polarization in Figs. (0), and 
(^. In Fig. (0), we explore the influence of ej on the relative power spectrum and the behavior is qualitatively similar 
to the Newtonian case. We explore, in Fig. (||), the effect of changing values for k by varying 7771,7712 and 77 after 
fixing the value of et- We see that the behavior is similar to the Newtonian case when we vary values of k for a given 
e. This is required as at 2PN order, for a given et, fc is uniquely determined by 7771,7772 and 77. However, there are 
quantitative differences in that positions and strengths of various harmonics are different in the Newtonian and the 
2PN cases. 

A quantitative comparison between the spectral analysis with Newtonian and 2PN motion is presented in Figs. IM) 
and (|I^. Note that we can perform this comparison as we are using scaled polarization waveforms, and H^. 
In Fig. (H), we plot both H^{1) and its power spectrum using Newtonian and 2PN accurate orbital motion. We 
choose the arbitrary parameter k, introduced in the Newtonian case to match the 2PN accurate k associated with 
the generalized quasi-Keplerian representation. In this manner, we force orbital elements for Newtonian and 2PN 
dynamics to be the same. Though, qualitatively similar, the plots for the Newtonian and the 2PN orbital motion are 
quantitatively different in that strengths of spectral lines are different by a few parts in thousand in most cases. 

Finally, in Fig. (0), we plot Hl{l) and its power spectrum as a function of orbital inclination angle for Newtonian 
and 2PN accurate orbital motion. The value of k for the Newtonian runs are again chosen so that it is comparable 
to the actual 2PN accurate k value. It is clear from these figures that inclusion of PN corrections to orbital motion 
changes distribution of spectral lines, though the position of the dominant (maximum amplitude) harmonic is roughly 
the same. This figure also shows how the orbital inclination angle i slowly modulates the spectral lines for Newtonian 
and 2PN orbital motion. 



IV. CONCLUSIONS 



A. Summary of results 



In this paper we have computed all the 'instantaneous' 2PN contributions to and hx for two compact objects 
of arbitrary mass ratio moving in elliptical orbits, using 2PN corrections to hfj'" and the generalized quasi-Keplerian 
representation for the 2PN motion. The expressions for and hx obtained here represent gravitational radiation 
from an elliptical binary during that stage of inspiral when orbital parameters are essentially the same over a few 
orbital periods, in other words when the gravitational radiation reaction is negligible. We investigate the effect of 
eccentricity, advance of periastron and orbital inclination on the power spectrum of the Newtonian part of and ■ 
The 2PN accurate generalized quasi-Keplerian representation is used in conjunction with two angular variables I and 
A chosen to facilitate the subsequent analysis of the waveform evolving under gravitational radiation reaction. These 
expressions thus form the first step in the direction of obtaining 'ready-to-use' theoretical templates for inspiralling 
compact bodies moving in quasi- elliptical orbits. 



B. Future directions 



There are several issues that remain open for further investigation. We list them below. 
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1. The next natural step is to obtain evolving hx and when lowest order radiation reaction effects are included 
in the evolution of orbital elements. This is currently under investigation 

2. There are tail contributions to and hx appearing at 1.5PN and 2PN orders. Though the formal expressions for 
tail terms are available in |^ , |lo| , they need to be written down in a form similar to 'instantaneous' contributions 
to and hx presented in this paper. 

3. After computing 'ready-to-use' search templates for inspiralling binaries in 'quasi-elliptical' orbits, one will 
be able to address a variety of data analysis issues related to the observations of gravitational radiation from 
eccentric binaries in great detail. These could include defining a 'restricted post-Newtonian' waveform, extending 
to 2PN accuracy the effect of eccentricity on detection discussed in [ p9p^ ] where currently the orbital dynamics 
is restricted to leading Newtonian order only. 

4. Finally, the analysis of the present paper may be useful for detecting continuous gravitational waves from known 
sources in binaries. Recent analysis [ p5[ employing a Keplerian representation for the binary's orbital motion 
indicates that the computational cost required to search such sources is affordable. The present analysis and 

may be crucial for including the relevant relativistic effects. 

Note added in proof: We observe that the results of our spectral analysis are in agreement with po| for low values 
of k. Since the effect of the periastron precession on the amplitude of spectral lines is not fully taken into consideration 
in [Q, the strength of the spectral lines in our analysis differs from theirs, for high values of k. It should be noted 
that for a given k, both methods give the same frequency shift for spectral lines, for all values of k. 
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'i Line of nodes 

FIG. 1. The orientation of unit vectors, which defines x and + polarization waveforms. The unit vectors p and q are 
the gravitational wave's principal axes with q = N x p. Note that N is a unit vector lying along the radial direction to the 

detector and p lies along the line of nodes. The Newtonian angular momentum vector L = ^ r x v is normal to the orbital 
plane and helps to define orbital inclination angle i. In this paper, the origin ioi 4> = \ + W is +ve x-axis, hence it is related 
to 4>Biww by = 0BIWW — f . 
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FIG. 2. Plots for scaled GW polarization waveform, ff" as a function of the mean anomaly, I and the corresponding 
normalized relative power spectrum using Newtonian orbital motion, for various values of eccentricity e. Note in H'^ (l) a 
'burst' of GW emission near periastron passage and a shift in the position of the dominant harmonic in the power spectrum as 
e increases. In the Fourier domain, the former results in a broad frequency rich peak. In all panels, the (arbitrary) periastron 
precession constant and the orbital inclination angle take values 0.1 and ^ respectively. 
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FIG. 3. The configuration is similar to Fig.^ but in tlie panels, k, the (arbitrary) periastron precession constant is varied 
for fixed eccentricity e — 0.5 and orbital inclination angle « = f. Note the splitting and shifting of spectral lines from integer 
multiple values of fr as k is increased. 
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FIG. 4. Plots of scaled GW polarization waveform, H^. as a function of mean anomaly, I and its corresponding normalized 
relative power spectrum for Newtonian orbital motion, when the orbital inclination angle i is varied. In all frames, eccentricity 
e = 0.5 and periastron precession constant, k = 0.1. 
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FIG. 5. The ratio of the total power measured iu x and + polarization for Newtonian motion as a function of orbital 
inclination angle i for various values of periastron precession constant k and eccentricity e. From the plots, it is clear that the 
ratio is independent of the orbital elements like e and k. 
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FIG. 6. Plots of A2c(0,A2s(0 as a function of mean anomaly, I and their relative power spectra constructed using Eq. 



(3.7) and directly using H'^{1). The orbital motion is Newtonian accurate and we employ certain specific values of A: = 0,0.5 
and 0.25. The relative power spectra, plotted in third column, are numerically identical, hence indistinguishable. Note that 
plots in second column are always 2n periodic, while those in the first column are 2 7r,2 7r and 47r periodic for k — 0, k — 0.5 
and k — 0.25 respectively 
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FIG. 7. Plots for scaled GW polarization waveform, (l) and corresponding normalized relative power spectrum using 
2PN accurate orbital motion, for various values of et- Unlike in Fig. the value of k cannot be independently chosen since it 
is uniquely determined by the values of mi,m2,n and et. The observations from plots are similar to Fig. as variation in k 



is small compared to that in et. In all frames, mi 
radians per second. 



m2 = 1.4 Mq, orbital inclination angle « = f- and mean motion n = 100 



37 



=0.4, i =Jt/3, (|)(,=7i/500} 




Mean anomaly, 1 




0.4 
0.2 


0.4 
0.2 


0.4 

0.2 — 













n= 100, k =0.142 

1 , . 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 











5 


10 1 




1 






n = 200, k =0.253 
1 , . 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



15 e 

O 

2; 



10 



15 



1 






n = 400, k =0.470 
III, 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 







5 


10 15 



Frequency in units of 



FIG. 8. The configuration is similar to Fig. |^ but in the panels we vary masses and mean motion to get the diflterent values 
of fe, rather than e. For all plots ej — 0.4 and i = §• Qualitatively, conclusions are similar to Fig. ^. Unit for n will be radians 
per second. 
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FIG. 9. Plots of scaled GW polarization waveform, Hl{l) and corresponding normalized relative power spectrum for 
Newtonian and 2PN accurate orbital motion. Given et, wo vary values for mi, m-z and n, to make k values the same in the 
2PN and the Newtonian accurate orbital motion. In all frames, orbital inclination angle is ^. Panels in the 2nd and 4th rows 
are for Newtonian orbital motion, whereas panels in the 1st and 3rd rows are for 2PN accurate orbital motion. We see minor 
quantitative differences in the position and strength of spectral lines at these two orders. 
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FIG. 10. Plots for the scaled GW polarization waveform, as a function of mean anomaly, I and the corresponding 
normalized relative power spectrum using 2PN and Newtonian accurate orbital motion. We vary the orbital inclination angle, 
i keeping other orbital elements constant at Newtonian and 2PN level. Conclusions are similar to Fig. ^. 
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